Origin of the Fermi Bubble 

K.-S. Cheng^ D. O. Chernyshov^'^'S, V. A. Dogiel^'^, C.-M. Ko^, W.-H. Ip^ 
^ Department of Physics, University of Hong Kong, Pokfulam Road, Hong Kong, China 
^ I.E.Tamm Theoretical Physics Division of P.N.Lebedev Institute of Physics, Leninskii pr. 

53, 119991 Moscow, Russia 
^ Institute of Astronomy, National Central University, Chung-Li 32054, Taiwan 



Received 



accepted 



- 2 - 
ABSTRACT 

Fermi has discovered two giant gamma-ray-emitting bubbles that extend 
nearly lOkpc in diameter north and south of the galactic center (GC). The ex- 
istence of the bubbles was first evidenced in X-rays detected by ROSAT and 
later WMAP detected an excess of radio signals at the location of the gamma- 
ray bubbles. We propose that periodic star capture processes by the galactic 
supermassive black hole, Sgr A*, with a capture rate 3 x lO^^yr"^ and energy 
release ~ 3 x lO^^erg per capture can produce very hot plasma ~ lOkeV with 
a wind velocity ~ lO^cm/s injected into the halo and heat up the halo gas to 
~ IkeV, which produces thermal X-rays. The periodic injection of hot plasma 
can produce shocks in the halo and accelerate electrons to ~TeV, which produce 
radio emission via synchrotron radiation, and gamma-rays via inverse Compton 
scattering with the relic and the galactic soft photons. 

Subject headings: Galaxy: halo - galaxies: jets - radiation mechanisms: non-thermal - 
black hole physics 
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Introduction 



Observations reveal many evidences of unusual processes occurring in the region of GC. 
For instance, the enigmatic 511 keV annihilation emission discovered by INTEGRAL (see 



e.g 



Knoedlseder et al 



20051 ) whose origin is still debated, the hot plasma with temperature 
about 10 keV which cannot be confined in the GC and, therefore, sourc es with a power 
about 10^^ erg s~^ are required to heat the plasma ( iKoyama et al.ll2007| ). In fact plasma 



outflows with the velo city > 100 km s ^ are observed in the nucle us regions of our Galaxy 



( ICrocker et al 



2010aJ) and Andromeda (iBogdan &: Gilfanov 



20101). Time variations of the 



6.4 keV line and X-ray continuum emission observed in the direction of molecular clouds 
in the GC which are suppose d to be a reflect i on of a giant X-ray flare occurred se veral 



hundred years ago in the GC (llnui et al. 



2009 



Ponti et al. 



2010 



Terrier et al. 



20101 ). HESS 



observations of the GC in the TeV energy range indicated an explosive injection of CR 
there which might be associated with the supermassive black hole Sgr A* (e.g. Aharonian et 
al. 2006). 



Recent analysis of Fermi LAT data (see 



Su et al. 



2010 



Dobler et al 



20101) discovered 



a new evidence of the GC activity. They found two giant features of gamma-ray emission 
in the range 1 -100 GeV, extending 50 degrees (~ 10 kpc) above and below the Galactic 
center - the Fermi Bubble (FB). They presented a list of mechanisms that may contribute 
to the energy release and particle production necessary to explain the gamma-ray emission 
from the bubble. They noticed, however, that most likely the Fermi bubble structure was 
created by some large episode of energy injection in the GC, such as a past accretion event 
onto the central supermassive black hole (SMBH) in the last ~ 10 Myr. They cast doubt 
on the idea that the Fermi bubble was generated by previous starburst activity in the GC 
because there was no evidence of massive supernova explosions (~ 10^ — 10^) in the past 
~ 10'' yr towards the GC. Besides, these supernova remnants should be traced by the line 
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emission 



of radio active ^^Al. Obser vations do not show significant concentration of ^^Al line 



towards the GC (jPiehl et al 



20061). 



Crocker fc Aharoniaru (120101 ) and 



used by 



Su et al. 



Crocker et al. 



fl2010a 



b|) argued that the procedure 



(120101 ) did not remove contributions of CR interaction with an ionised gas, 
then gamm a-rays could be produc e d by protons interaction with the fully ionised plasma 



in the halo. ICrocker fc AharonianI ( l2010l ) also argued that the lifetime of these protons can 
be very long because the plasma is extremely turbulent in this region therefore protons 
could be trapped there for a time scale Tpp > 10^° yr, then the observed gamma-rays can be 
explained with the injected power of SN ~ 10^^ erg s~^. 



In this letter we propose t hat the FB emission 



which have been developed by 



Cheng et al. 



(12006 



may 



resul t from star capture p ro cesses , 



20071 ) and 



Dogiel et al. 



(l2009a 



d|) to 



explain a wide range of X-ray and gamma-ray emission phenomena from the GC. 



2. Observations 



The procedure of sep aration o f the b ubble emission from the total diffuse emission of 



the Galaxy is described in 



Su et al. 



3. 



(120101 ). It is important to note that the bubble structure 



is seen when components of gamma-ray emission produced by cosmic ray (CR) interaction 
with the background gas, i.e. by CR protons (vr" decay) and electrons (bremsstrahlung) are 



removed. 



Su et al 



(|2010[ ) concluded the the bubble emission was of the inverse Compton 
(IC) origin generated by relativistic electrons. Here we summarize the multi-wavelength 
observational constraints of the FB: 



The spectral shape and intensity of gamma-rays are almost constant over the bubble 
region that suggests a uniform production of gamma-rays in the FB. The total 
gamma-ray flux from the bubble at energies > 1 GeV is ~ 4 x 10'^'' erg s"^ 
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and the pho ton spectrum o f gamma-rays is power-law, (iA^^/ cii?^ oc E^'^ for the range 



1-100 GeV flSu et al. 



(I2OI0I)); 



In the radio range the bubble is seen from the tens GHz WMAP data as a microwave 
residual spherica 



( iFinkbeiner 



excess ("the microwave haze") above the GC ~ 4 kpc in radius 
2OO4J ) ■ Its power spectrum in the frequency range 23 - 33 GHz is described 
as power-law, $1, oc u"^-^. For the magnetic field strength ~ 10 fiG the energy 
range of electrons responsible for emitting these radio waves is within the range 20 - 
30 GeV and their spectrum is dN^/dEf. oc E^ 



'2. 

■'e ) 



• The ROS AT 1.5 keV X-ray data clearly s howed the characteristic of a bipolar 



structure (IBland-Hawthorn fc Cohen 



2OO3I ) that aligned well with the edges of the 



Fermi bubble. The ROSAT structure is explained as due to a fast wind which drove a 
shock into the halo gas with the velocity Vsh ~ lO^cm s~^. This phenomenon requires 
an energy release of about 10^^ ergs at the GC and this activity should be periodic on 
a timescale of 10 — 15 Myr. 

The similarities of the morphology of the radio. X-ray and gamma-ray structures 
strongly suggest their common origin. 



Su et al. 



( I2OIOI ) gamma-rays are produced 



In the case of electron (leptonic) model of 
by scattering of relativistic electrons on background soft photons, i.e. relic, IR and optical 
photons from the disk. 



3. The bubble origin in model of a multiple star capture by the central black 

hole 

In this section we present our ideas about the origin of the Fermi Bubble in the 
framework of star capture by the central SMBH. The process of gamma-ray emission from 
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the bubble is determined by a number of stages of energy transformation. Each of these 
stages actually involves complicated physical processes. The exact details of these processes 
are still not understood very well. Nevertheless, their qualitative features do not depend 
on these details. In the following, we only briefly describe these processes and give their 
qualitative interpretations . We begin to descr ibe processes of star capture by the central 



black hole as presented in iDogiel et al.l ( ]2009dl ). 



3.1. Star capture by the central black hole 



As observations show, there is a supermassive black hole (Sgr A*) in the center of 
our Galaxy with a mass of ~ 4 x 10^ Mq . A total tidal disruption of a star occurs when 
the penetration parameter h^^ ^ 1, where h is the ratio of the periapse distance Tp to the 
tidal radius Rt- The tidal disruption rate Vg can b e approxiraated to within an order of 



magnitude i/^ ~ 10 ^ — 10 ^ yr ^ (see the review of 



Alexander 



20051) 



The energy budget of a tidal disruption event follows from analysis of star matter 
dynamics. Initially about 50% of the stellar mass becomes tightly bound to the black hole , 
while the remainder 50% of the stellar mass is forcefully ejected (see, e.g. I Aval et al.l 120001 ). 
The kinetic energy carried by the ejected debris is a function of the penetration parameter 
h^^ and can significantly e xceed that releas ed by a normal supernova (~ 10^^ erg) if the 



orbit is highly penetrating (lAlexander 



20051 ) 



ly ~ 4 X 10^2 



Mo 



Rq 



106 ) Vo.i 



erg . 



(1) 



Thus, the mean kinetic energy per escaping nucleon is estimated as i^csc ~ 
42 [fiY^ (^) (j^^^f^y^ (o^)~^ MeV, where r]M, is the mass of escap- 

ing material. From W and Vg we obtain that the average power in the form of a flux of 
subrelativistic protons. If ~ 3 x 10^^ erg and z/^ ~ 3 x 10"^ yr~^, we get ~ 3 x 10"^° 
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erg s . 

In iDogiel et al.l (l2009d ) we described the injection spectrum of protons generated 
by processes of star capture as monoenergetic. This is a simphfication of the injection 
process because a stream of ch arged particles should be influenced by different plasma 



instabilities, as it was shown by 



Ginzburg et al. 



(j2004| ) for the case of relativistic jets. At 



first the jet material is moved by inertia. Then due to the excitation of plasma instabilties 
in the flux, the particle distribution functions, which were initially delta functions both in 
angle and in energy, transforms into complex angular and energy dependencies. 



3.2. Plasma heating by subrelativistic protons 



dEp 
dt 



Subrelativistic protons lose their energy mainly by Coulomb collisions, i.e. 

= — In A, where Vp is the proton velocity, and In A is the Coulomb logarithm. 



In this way the protons transfer almo st all their e n ergy t o the background plasma and 



heat it. This process was analysed in iDogiel et al.l (i2009cu dl). For the GC parameters the 
average time of Coulomb losses for several tens MeV protons is several million years. The 
radius of plasma heated by the protons is estimated f rom the diffusion eq uation describing 



2009bh . This radius is 



propagation and energy losses of protons in the GC fiDogiel et al. 
about 100 pc. The temperature of heated plasma is determined by the energy which these 
protons transfer t o the background ga s. For W ~ 10^° — lO^^erg s~^ the plasma temperature 



is about 10 keV ( iKoyama et al. 



20071 ) just as observed by Suzaku for the GC. The plasma 
with such a high temperature cannot be confined at the GC and therefore it expands into 
the surrounding medium. 
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e.g 



3.3. The hydrodynamic expansion stage 

Hydrodynamics of gas expansio n was described in many monographs and reviews (see 



Bisnovatyi-Kogan Sz Silich 



19951 ). As the time of star capture may be smaller than the 
time of proton energy losses, we have almost stationary energy release in the central region 
with a power W ^ 3 x 10*^'' erg s~^. This situation is very similar to that described by 



Weaver et al. 



(119771 ) for a stellar wind expanding into a uniform density medium. The 
model describes a star at time t=0 begins to blow a spherically symmetric wind with a 
velocity of stellar wind V^, mass-loss rate dM^/dt = M^, and a luminosity = 
which is analogous to the power W produced by star capture processes. Most of the time 
of the evolution is occupied by the so-called snowplow phase when a thin shock front is 
propagating through the medium . The shock is expanding as 

Rskit) = a i^^j (2) 

where po = nonip and a is a constant of order of unity. The velocity distribution inside the 
expanding region u{z) is nonuniform. 

Our extrapolation of this hydrodynamic solution onto the Fermi bubble is, of course, 
rather rough. First, the gas distribution in the halo is nonuniform. Second, the analysis 
does not take into account particle acceleration by the shock. A significant fraction of the 
shock energy is spent on acceleration that modify the shock structure. Nevertheless, this 
model presents a qualitative picture of a shock in the halo. 



3.4. Shock wave acceleration phase and non-thermal emission 

The theory of particle acceleration by shock is described in many publications. This 
theory has been developed, and bulky numerical calculations have beeen performed to 
calculate spectra of particles accelerated by supernova (SN) shocks and emission produced 
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by accelerated electrons and protons in the range from radio up to TeV gamma-rays (see 



e.g. 



Berezhko fc Voelk 



20101). Nevertheless many aspects of these processes are still unclear. 



For instance, the ratio of el ectrons to protons accelerated by shocks is not reliably estimated 



see 



Bykov fc Uvarov 



I99I . 



We notice that the energy of shock front expected in the GC is nearly two orders of 
magnitude larger than that of the energy released in a SN explosion. Therefore, process of 
particle acceleration in terms of sizes of envelope, number of accelerated particles, maximum 
energy of accelerated particles, etc. may differ significantly from those obtained for SNs. 
Below we present simple estimations of electron acceleration by shocks. 



2 

e ) 



The injection spectrum of electrons accelerated in shocks is power-law, Q{Ee) oc 
and the maximum energy of accelerated electrons can be estima ted from a kinetic equation 



19901), E^^a. ^ vlJDP, where 



describing their spectrum at the shock (iBerezinskii et al. 
Vsh ~ lO'^cm/s is the velocity of shock front, D is the diffusion coefficient at the shock front 
and the energy losses of electrons (synchrotron and IC) are: dEe/dt = —I3E^. Recall /3 is 
a function of the magnetic and background radiation energy densities, /3 ~ waTc/ {mf.c^Y' , 
where ctt is Thompson cross section and w = Wph + wb is the combined energy density of 
background photons Wph and the magnetic energy density wb respectively. It is difficult to 
estimate the diffusion coefficient near the shock. For qualitative estimation, we can use the 
Bohm diffusion (~ rL{Ee)c), where is the Larmor radius of electrons. Using the typical 
values of these parameters, we obtain E^ax ~ 1 TeV VsB]/^wZli^ , where Vg is the shock 
velocity in units of 10®cm/s, is the magnetic field in the shock in units of 10~^G and 
w_i2 is the energy density in units of 10~^^erg/cm^. 

The spectrum of electrons in the bubble is modified by processes of energy losses and 

escape. It can be derived from the kinetic equation 

d fdE\ N ^^^^ 
lE [HV + T = «<^) 
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where dE/dt = PE'^ + Vu{z)E describes the inverse Compton, synchrotron and adiabatic 
(because of wind velocity variations) energy losses, T is the time of particle escape from 
the bubble, and Q{E) = KE^'^9{Emax — E) describes particles injection spectrum in the 
bubble. As one can see, in general case the spe ctrum of electrons in the bubble cannot be 



Suetal 



fcOlOl ). The spectrum of electrons 



described by a single power-law as assumed by 
has a break at the energy Ef, ~ l/PT where T is the characteristic time of either particle 
escape from the bubble or of the adiabatic losses, e.g. for Vu = a the break position follows 
from T ~ l/a. By solving equation 3, we can see that the electron spectrum cannot be 
described by a single power-law even in case of power-law injection (see eg. Berezinskii et 
al., 1990). 

The distribution of background photons can be derived from GALPROP program. 
The average energy density of background photons in the halo are Wo = 2 eV cm~^ for 
optical and wjr = 0.34 eV cm~^ for JR. These background soft photon energy densities are 
obviously not negligible in comparing with wcmb = 0.25 eV cm~^ for the relic photons 
and also comparable with the magnetic energy density (~ l{H/5 x lO^^GyeV/cm^). The 
expected IC energy flux of gamma-rays and synchrotron radiation emitted from the same 
population of electrons described above are sh own in Fig. [J for different values of Eb 



and Emax- The Klein-Nishina IC cross-section ( iBlumenthal fc Gould 



mm is used. The 



observed spectrum o: 



Dobler fc Finkbeinerl fl2008[ ) and 



radio emiss ion in the range 5-200GHz and gamma-rays are taken from 



Su et al. 



(I2OIOI ) respectively. The inverse Compton gamma 
ray spectrum is formed by scattering on three different components of the background 
photons. When these three components are combined (see Fi g. Hb), they mi mic a photon 



spectrum E ^ and describe well the data shown in Fig. 23 of 



Su et al. 



mm . We want to 



remark that although a single power law with the spectral indexes in between 1.8 and 2.4 in 
the energy range of electrons from 0.1 to 1000 GeV can also explain both the Fermi data as 
well as the radio data as suggested by Su et al. (2010). Theoretically a more complicated 



electron spectrum will be developed when the cooling time scale is comparable with the 
escape time even electrons are injected with a single power law as shown in equation 3. 



3.5. The thermal emission from heated plasma 

In our model there is lOkeV hot plasma with power W ^ 3 x lO'^^erg/s injected 
into bubbles. Part of these energies are used to accelerate the charged particles in the 
shock but a good fraction of energy will be used to heat up the gas in the halo due to 
Coulomb collisions. The temperature of halo gas can be estimate as nd^kT Wd/v^j 
which gives kT ~ 1.5(1^/3 X 10^^erg/s)v^^{d/5kpc)-^{n/10-^)keV. The thermal radiation 
power from the heated halo gas is simply given by Lth = 1.4 x 10~'^'^ rieriiZ'^T^^'^er g/cm'^ 



(iRybicki &: Lightman 



19791). By using kT=1.5keV, Ug = rii = 10 ^cm ^ and Z=l, we find 



~ lO^^erg/s. 



4. Discussion 



The observed giant structure of FB is difficult to be explained by other processes. We 
suggest that periodic star capture processes by the central SMBH can inject ~ 3 x lO'^^erg/s 
hot plasma into the galactic halo. The hot gas can expand hydrodynamically and form 
shock to accelerate electrons to relativistic speed. Synchrotron radiation and inverse 
Compton scattering with the background soft photons produce the observed radio and 
gamma-rays respectively. Acceleration of protons by the same shock may contribute to 
Ultra-high energy cosmic rays, which will be considered in future works. 

It is interesting to point out that the mean free path of TeV electrons A ~ \/D/ /3Ef, ~ 
bOD^l^T^^'^pc, where D2S and are the diffusion coefficient and cooling time for TeV 
electrons in units of lO^^cm^/ s and lO^yrs respectively. This estimated mean free path is 
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much shorter than the size of the bubble. In our model the capture time is once every 
~ 3 X lO^yr, we expect that there is about nearly 100 captures in 3 million years. Each of 
these captures can produce an individual shock front, therefore the gamma-ray radiation 
can be emitted uniformly over the entire bubble. 

Furthermore we can estimate the shape of the bubble, if we simplify the geometry of 
our model as follows. After each capture a disk-like hot gas will be ejected from the GC. 
Since the gas pressure in the halo {n{r)kT ~ 10-^''(n/3 x 10-^cm-=^)(T/3 x 10^ir)erg/cm3) 
is low and decreases quadratically for distance larger than 6kpc(Paczynski 1990), wc can 
assume that the hot gas can escape freely vertically, which is defined as the z-direction and 
hence the z-component of the wind velocity = constant or z — Vy^^t- The ejected disk 
has a thickness Az = v^^tcap, where tcap = 3 x lO^yr is the capture time scale. On the 
other hand the hot gas can also expand laterally and its radius along the direction of the 
galactic disk is given by x{t) — v^xt + a^o ~ ''^wxt, where Xq ~ 100pc(cf. section 3.2). When 
the expansion speed is supersonic then shock front can be formed at the edge of the ejected 
disk. In the vertical co-moving frame of the ejected disk the energy of the disk is AE, 
which is approximately constant if the radiation loss is small. The energy conservation gives 
AE — \mv^^ with m — mo + ms{t) — mo + nx'^Azp, where mo ~ lAEjv^ is the initial 
mass in the ejected disk, m^ is the swept-up mass from the surrounding gas when the disk 
is expanding laterally, and p = nipn is the density of the medium surrounding the bubble. 
Combing the above equations, we can obtain AE — |[mo + Tr{vwxty Azplv^^^. There are 
two characteristic stages, i.e. free expansion stage, in which v^x ~constant for mo > ms(i) 
and deceleration stage for mo < ^^{t). The time scale switching from free expansion to 
deceleration is given by mo = rns{ts) or tg = ^/mo/nAzpvl^^. In the free expansion stage, we 
obtain x — ^z ~ 2; for x < Vyjxts — Xg. In the deceleration stage, AE ^ ^7r{vyjxty^zp'^wx^ 
we obtain {xjv^t,) = (^^i^^j (^j ~ 0.9 j , we have approximated 
Vwx ~ Vyuz ^ Vw ^ lO^cm/s, AE = 3 x lO^^ergs and p/nip = 3 x 10~^cm~"^. The switching 
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from a linear relation to the quadratic relation takes place at Zg ~ v^tg ~ 300pc. The 
quasi-periodic injection of disks into the halo can form a sharp edge, where shock fronts 
result from the laterally expanding disks with quadratic shape, i.e. z ^ x^. In fitting the 
gamma-ray spectrum it gives E}, ~ 50GeV, which corresponds to a characteristic time scale 
of either adiabatic loss or particle escape ~ ISMyrs. By using equation 2, the characteristic 
radius of FB is about 5kpc, which is quite close to the observed size of FB. 

We thank the anonymous referee for very useful suggestions and Y.W. Yu for useful 
discussion. KSC is supported by a grant under HKU 7011/ lOp. VAD and DOC are partly 
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CMK is supported in part by the National Science Council, Taiwan, under grants NSC 
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1. — (a)Energy fluxes produced by Synchrotron and IC with £'^a^=lTeV and different 
(b)The gamma-ray spectra with Eh=50GeV and different Ej^ax- ■ 
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